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To shed light into the pairing mechanism of possible spin-triplet superconductors (TMTSF^X 
and Sr2Ru04, we study the competition among various spin singlet and triplet pairing channels in 
the Hubbard model by calculating the pairing interaction vertex using the ground state quantum 
Monte Carlo technique. We model (TMTSF^X by a quarter-filled quasi-one dimensional (quasi-f D) 
Hubbard model, and the 7 band of S^RuCU by a two dimensional (2D) Hubbard model with a band 
filling of ~ 4/3. For the quasi-lD system, we find that triplet /-wave pairing not only dominates over 
triplet p-wave in agreement with the spin fluctuation theory, but also looks unexpectedly competitive 
against d-wave. For the 2D system, although the results suggest presence of attractive interaction 
in the triplet pairing channels, the d-wave pairing interaction is found to be larger than those of the 
triplet channels. 
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I. INTRODUCTION 



Possible occurrence of spin-triplet pairing in super- 
conductors such as a ruthenate S^RuO^i* 2 * 3 ^ and or- 
ganic materials (TMTSF) 2 X (X=C10 4) PF 6 .„)&& has 
attracted much attention lately. A probable mechanism 
of spin-triplet pairing is due to ferromagnetic spin fluc- 
tuations, as in 3 He, as proposed in refQ for S^RuC^. 
In both of these materials, however, the mechanism of 
triplet pairing has been controversial because the spin 
fluctuation has turned out to be more like antiferromag- 
netic than ferromagnetic. Namely in (TMTSF^X, the 
nesting of the quasi-lD Fermi surface give rise to a 2fcp 
spin density wave phase lying right next to the super- 
conducting phase. In Sr2Ru04, spin fluctuations with a 
wave number (2ir/3, 2tt/3) arises due to the nesting of 
the quasi-lD a and (3 bands£ 

For (TMTSF) 2 X, the possibility of spin triplet p-wave 
pairing has been considered from the early dayS)2ii2iii 
but if we consider the 2kp spin fluctuations to be the 
origin of the pairing interaction, spin singlet d-wave- 
like pairing is more likely to occur as proposed by sev- 
eral authors piSiiiiii^ Recently, two of the present au- 
thors have proposed, based on a phenomelogical spin- 
charge fluctuation theory, that a combination of quasi- 
one-dimensionality, the coexistence of 2kp spin and 
charge fluctuations, and anisotropy in the spin fluctu- 
ations may lead to spin-triplet /-wave pairing, dominat- 
ing over d- and p-wave pairings^ A similar proposal has 
been made by Fuseya et alii In fact, the coexistence 
of 2k p spin and charge density wave has been reported 
experimentally." 1 ^* 1 ^ However, a recent renormalization 
group study by Fuseya et alm& has suggested that p-w&ve 
pairing dominates over /-wave pairing in quasi- ID sys- 



tems. 

For Sr2RuC>4 on the other hand, there has been a 
debate as to which one of the three bands plays the 
main role in the occurrence of superconductivity, namely 
the 2D 7 bands or the quasi-lD a-j3 bands. Several 
microscopic theories have focused on the antiferromag- 
netic spin fluctuations due to the a-(3 band Fermi surface 
nestingi 21 i 22 i 23 i 24 i 25 In refs l2lll23l in particular, the exper- 
imentally observed anisotropy^ in the spin fluctuations 
is the key for spin-triplet p-wave pairing to take place. 

On the other hand, Nomura and Yamada focused 
mainly on the 2D 7 bandii Applying third order per- 
turbation theory to the 2D Hubbard model with appro- 
priate band structure and band filling, they showed that 
triplet pairing superconductivity takes place, dominating 
over singlet pairing. A multiband calculation has recently 
been performed along this liner 28 Also, a renormalization 
group study has been performed for the Hubbard model 
having nearest neighbor hopping, where spin triplet pair- 
ing is shown to be the leading instability near but away 
from the van Hove filling iSi By contrast, if we apply ran- 
dom phase approximation (RPA) type theories like the 
fluctuation exchange method (FLEX)SI to the 2D Hub- 
bard model, triplet pairing not only does not occur in a 
realistic temperature regime,— but also does not even 
dominate over singlet pairing at the filling corresponding 
to the 7 band of Sr 2 Ru04 (see sectio JV"D|) . 

The purpose of the present study is to shed light 
into these controversies. We model (TMTSF^X and 
Si'2Ru04 by quasi-lD and 2D Hubbard models, respec- 
tively. Using the ground state quantum Monte Carlo 
technique, we calculate the pairing interaction vertices 
for possible pairing symmetries. 
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(see Fig0. The pairing symmetry is determined so as to 
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FIG. 1: p, d, and /-wave gap functions along with the Fermi 
surface of the quasi ID system and the nesting vector Q. The 
arrows denote a pair scattering process having a momentum 
transfer Q. 
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FIG. 2: A generic phase diagram for quasi-lD systems. 



II. SPIN-CHARGE FLUCTUATION THEORY 

Before going into the QMC calculation, here we sum- 
marize the spin-charge fluctuation theory of quasi-lD sys- 
tems given in refll6t The pairing interactions mediated 
by spin and charge fluctuations are generally given as 



Vs(q) = \v^ z {q) + V+-(q) l -V ch (q) 



for singlet pairing, 
V t ±(q) 



^V z p z (q)--V ch (q) 



(1) 



(2) 



for triplet pairing with total S z = ±1 (d- vector d _L z), 
and 

V t \\ (q) = \v£{q) V+-(q) \v ch {q) (3) 



positive and large, where V a (k — k') is one of the pairing 
interactions iJTJ-GJ and <Xfc) is the gap function. When 
V a (q) peaks around the nesting vector q — Q, the con- 
tribution around k — k' — Q becomes dominant in the 
summation of eqn.Q), so that the sign of the quantity 
V a (Q)<j)(k)(f){k + Q) has to be negative in order to have 
superconductivity. 

Let us now consider p, d, /-wave gap functions shown 
in Fig^ along with the Fermi surface. Although the 
d- (/-) wave gap is not d (/) in the strict sense of the 
word, here we use this terminology in the sense that the 

gap changes sign asH 1 — (H 1 1 — ) along 

the Fermi surface. Since the ci-wave gap function sat- 
isfies <j){k)4>{k + Q) < 0, a positive pairing interaction 
V S (Q) is necessary for superconductivity to occur in this 
channel. Similarly V t ±(Q) < or V t \\{Q) < is nec- 
essary for /-wave pairing. If we assume V s z p z > , 
|Ft_L(Q)| > |Vf||(Q)|, so that /-wave is more likely to oc- 
cur in the d _L z channel. Since the number of nodes 
of the gap intersecting the Fermi surface is the same 
between d and / due to the quasi-one-dimcnsionality 
the competition between the two is almost solely dom- 
inated by the absolute value of the pairing interaction. 
Thus the boundary between d and /-waves is roughly 
given by Vt±(q) — —V s (q), which corresponds to the line 
V c h = ■ I n particular, when the system has spin 
rotational symmetry, the two symmetries become nearly 
degenerate when spin and charge fluctuations have about 
the same strength^ 

On the other hand, the p-wave gap function satisfies 
<j)(k)(j)(k + Q) < 0, so that a positive triplet pairing in- 
teraction is necessary. This can be the case only for 
d || z and only when V c h/V s Z p and V+~ /V s z p z are both 
small iSiiSSiSi Since the p-wave gap function is nodeless 
on the Fermi surface as opposed to the case of d- and 
/-waves, this will be the dominant pairing if the pairing 
interaction is positive and strong enough. 

From the above phenomelogical argument, we obtain 
a generic phase diagram given in Fig|3 



for triplet pairing with S z = (d \\ z). Here V z p z , V s ^~ , 
and V c h (all positive) are contributions from the longi- 
tudinal spin, transverse spin, and charge fluctuations, 
respectively. We assume that the fluctuations are en- 
hanced at q = Q, where Q is the nesting vector of the 
Fermi surface. In the case of (TMTSF) 2 X, Q bridges the 
two disconnected pieces of the quasi ID Fermi surface 



III. FORMULATION OF THE QMC STUDY 



In this section, we formulate the ground state QMC 
study. 
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FIG. 3: The on-site Hubbard model for 
system, and (b) the 2D system. 



'a) the quasi-lD 



A. On-site Hubbard model 

As a model for (TMTSF) 2 X, we consider the Hubbard 
model 



F=-5>i(4 



on an N x x N y quasi-lD lattice shown in Fig|3fa). We 
set t x <C t y , where t x (taken as unity throughout the 
study) and t y are hopping integrals in the x and y di- 
rections, respectively. Since the energy band is much 
more dispersive in the x direction than in the y direc- 
tion, we take more sites in the former than in the latter 
(N x ^S> N y ). As for the electron-electron interaction, we 
assume only the on-site repulsion U for simplicity, and 
refer to the model as the on-site Hubbard model. The 
reason for this is to make the QMC calculation feasi- 
ble. Note that, U enhances only the spin fluctuation 
(Xp 2 = Kp" ^ ^ch)j so that d-wave pairing is expected 
from the spin fluctuation theorjiiSii 3 . as seen from the 
phase diagram. We take U = 2 throughout the study. 
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FIG. 4: The pairing channels considered for (a) the quasi- 
lD system, and (b) the 2D system. The signs denote g a (S), 
where S is measured from the center site. 



The band filling n = (number of electrons N e /number of 
sites TV) is taken close to quarter filling, i.e., n ~ 0.5. 

As for Sr2RuC>4, we focus on the 7 band and model 
it by the 2D single band Hubbard model having nearest 
and next nearest neighbor hoppings, t and t' respectively, 
as shown in Fig|3Jb). t' has been estimated to be 0.3-0.4 
in units of |t|. 33 i 34 i 35 Here we take t' = 0.4. The band 
filling is taken to be close to n = 4/3. 



B. Calculation conditions 

We adopt the auxiliary field quantum Monte Carlo 
method for the ground state^2i2S^2ii£ In this formalism, 
physical quantities are evaluated using the ground state 
wave function obtained as 

exp(-r J ff)|0 T ), 

where H is the Hamiltonian, r is a sufficiently large num- 
ber, and \4>t) is a trial state, r is sliced into a product of 
many At using Trotter-Suzuki decomposition. 36 In the 
actual calculation, we have taken (3 ~ 10, and At < 0.1, 
which turns out to be sufficient for a relatively small value 
of U = 2. The correlation functions are evaluated as the 
average of the calculation results of about 10 independent 
Monte Carlo runs, each of which consists of few thousand 
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MC sweeps. The error bars are estimated as the variation 
of the independent runs. 



C. Pairing interaction vertex 

We calculate the pairing interaction vertex for possible 
pairing channels, 39 which is defined as 

AP a (i,j) = ^2g a (S)g a (S') x 

6,5' 

(caCi + 5Tct T c} +5 , ; ) - 



± 



± 



\ c it c ]-|-) ( c i+<51 C j+<5'j, 

){ci+sic] l ) 

Here, the signsi takes + for singlet pairing, and — for 
triplet. g a (S) gives sign that depends on the pairing chan- 
nel a as given below. A positive interaction vertex in a 
certain pairing channel implies that there is an effective 
attractive pairing interaction in that channel between the 
renormalized quasi-particles. 

For the quasi-lD system, we consider the pairing chan- 
nels pi, p2, d and / as shown in Fig0Ja) along with 
the signs g a (S). It can be easily shown that the cor- 
responding gap functions in momentum space are pro- 
portional to sin(2fc x ), sin(fca;), cos(2fc x ), and sin(4fcj,), 
respectively. For the 2D system, three triplet pairings 
tpl, tp2, tp3 and a singlet pairing spl are considered as 
shown in FigQfb). The corresponding gap functions in 
momentum space are proportional to sin(k x ) (or equiv- 
alently sin(fc y )), sin(fc x + k y ) (or sin(fca; — k y )), sin(2fc x ) 
(or sin(2fc a )), and cos(fc x ) — cos(fep), respectively, spl is 
the so-called d x i_ y i pairing. 

Let us denote the i-th site as i = (i x ,i y ). For the 
quasi- ID system, we calculate 

AP a (r) = AP a (i x ,iy,j x ,i y ) 

ly 

(pairing correlation within each chain) as a function of 
distance in the x direction r = \i x — j x \. For the 2D 
system, we calculate 

as a function of the "distance" r. In both systems, we 
focus on the long range part AP Q (r) because supercon- 
ductivity is characterized by the long ranged behavior of 
the pairing correlation function. 



D. Open shell and closed shell fillings 

Since the energy scale of superconductivity in re- 
pulsively interacting systems is small, pairing correla- 



tions are known to be strongly affected by the en- 
ergy level spacing around the Fermi level in finite size 
systemsi 41 i 42 ' 43 Namely, since superconductivity occurs 
due to pair scattering between states near the Fermi level, 
the enhancement of the pairing correlation can be wiped 
out when the energies of the highest occupied levels and 
the lowest unoccupied levels are well separated. 

In order to cope with this problem, here we consider 
the following two situations, (i) Open shell fillings : the 
number of electrons are chosen so that the fc-points at 
the Fermi level is half-filled for the U = ground state. 
The highest occupied and the lowest unoccupied states 
have the same energy in this case. We tune the value of 
the hopping integrals (t y for the quasi-lD system and t' 
for the 2D system) appropriately so as to distribute the 
fc-points having the same energy as uniformly as possible 
along the Fermi surface. 

In the case of open shell fillings, the trial state is chosen 
as follows. The U = single particle states (exp(ifc • r) ) 
below the Fermi energy is completely filled with electrons. 
At the Fermi energy, the states of the form exp(ife • r) + 
exp(— ik ■ r) are filled only by up spin electrons, while 
the states of the form exp(ik ■ r) — exp(— ik ■ r) arc filled 
only by down spin electrons. We will use gray circles to 
denote the half-filled fc-points in the trial state, as shown 
in the inset of FigsEO |7| El andlTUTah 

(ii) Closed shell fillings : the number of electrons is 
chosen so that the Fermi energy for the U = ground 
state comes in between slightly separated energy levels. 
For the quasi-lD system, t y is tuned so that the energy 
difference between the highest occupied and the lowest 
unoccupied levels is within 0.01t x for U = 0. As for 
the 2D system, t y is taken to be slightly different from 
t x thereby lifting the degeneracy, by less than 0.01^, 
between the highest occupied levels k = (fei, fe) and the 
lowest unoccupied levels k — (k^, fci). 

In this case we take the U = ground state, which is 
a state with total spin Stot = 0, as the QMC trial state. 
Closed (open) circles will be used to denote the occupied 
(unoccupied) /c-points that lie within O.Olia; to the Fermi 
energy in the trial state, as shown in the inset of FigsEJ 
El and lTDf b). 

For open shell fillings, triplet pairing is expected to be 
favored because the electrons at the half-filled Fermi level 
tends to be in a high spin configuration. For closed shell 
fillings on the other hand, spin singlet pairing should be 
favored since we restrict ourselves to Stot = states by 
choosing the U = ground state as the trial state. This 
problem does make the study of singlet-triplet competi- 
tion less conclusive, but our results still give significant 
information concerning the d- vs. /-wave competition in 
the on-site Hubbard model as we shall see below. As for 
the competition within triplet pairings (p vs. /) on the 
other hand, we can safely give a robust conclusion. 
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IV. RESULTS AND DISCUSSIONS FOR THE 
QUASI ID SYSTEM 

A. Open shell filling 

We now show the results for the quasi ID system. We 
first present results for open shell fillings. The pairing in- 
teraction vertices AP S are shown for three different sets of 
parameters in Fig0 In all cases, APf is found to be posi- 
tive at large distances, while AP p i and AP P 2 are negative 
at most of the distances. This means that the spin-triplet 
pairing interaction is attractive in the /-wave channel 
and repulsive in the p-wave channels. This result is in 
contrast with the recent renormalization group study 20 , 
and is at least qualitatively consistent with the spin- 
charge fluctuation theory described in sectiorlTTl where 
V s f = V+- » V ch and thus V n (Q) = V t± {Q) < hold 
for the on-site Hubbard model. 

The interaction vertex for ti-wave pairing APd takes 
both positive and negative values, and we find APf > 
APd at most of the large distances^ This is somewhat 
surprising because when the system has spin rotational 
symmetry, /-wave can be competitive against <i-wave 
only when 2kp spin and 2kp charge fluctuations have 
about the same strength as far as the spin-charge fluctu- 
ation theory is concerned^ which is not the case for the 
on-site Hubbard models 

As mentioned before, open shell filling is expected to 
favor spin triplet pairing, so this result may be regarded 
as an artifact due to the open shell filling. However, as 
we shall see later, open shell filling does not necessarily 
result in AP t ripict > AP S j ng i c t. Another possible reason 
for / dominating over d is that effects beyond RPA-type 
theories may be important in quasi-lD systems. This is 
highly probable because Fermi liquid picture is known to 
be invalid at least in purely ID systems 



B. Closed shell filling 

Calculation results for a closed shell filling is shown 
in FigEl In this case, we have positive APd, smaller 
but still positive (or non-negative) APf for all distances, 
and clearly negative AP p . Thus, we may conclude that 
APf > AP p holds for both open and closed shell fillings. 
On the other hand, although the result of APd > APf 
apparently seems to be consistent with the spin fluctua- 
tion theory, this result may not be taken directly is since 
closed shell fillings selectively favors spin singlet pairing 
as mentioned above. We shall come back to this point 
later. 



xiQ- 3 3 
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FIG. 5: Pairing interaction vertex AP a plotted against dis- 
tance r for the quasi-lD system with open shell fillings. Pa- 
rameter values are taken as (TV, N e , n, t y )= (a) (28 x 6, 
84, 0.50, 0.158), (b) (32 x 6, 96, 0.50, 0.246) , (c) (24 x 8, 
100,0.52, 0.259) . The inset shows the half-filled fc-points just 
at the Fermi energy in the trial state of QMC, along with the 
Fermi surface of the corresponding infinite size systems. 



V. RESULTS AND DISCUSSIONS FOR THE 2D 
SYSTEM 

A. Open shell fillings 

Let us now move on to the 2D system with n ~ 4/3. 
We first look into the case with open shell fillings. In 
FigEIa), we focus on the interaction vertices for triplet 
pairings, where we find positive interaction vertices for 
all three channels. However, if we compare singlet and 
triplet channels in Fig[7fb), AP sp i > AP tp holds for all 
three triplet channels. A similar result is obtained also 
for a 14 x 14 system as seen in Fig|HJ Since the pair- 
ing forms are not optimized and there is a possibility 
of stronger triplet pairings between more separated sites 
than we have considered, we cannot conclude that sin- 
glet dominates over triplet. However, considering the fact 
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, (a) 



FIG. 6: A plot similar to Fig|^| but for a closed shell Pa- 
rameter values are taken as (N, N e , n, t y ) = (32 x 6, 98, 0.51, 
0.240). The inset shows the fc-points that lie within 0.01t x to 
the Fermi energy in the trial state of the QMC calculation, 
along with the Fermi surface of the corresponding infinite size 
systems. 



that this is a result for an open shell filling, we may say 
that d x 2_ y 2-wave is still a strong candidate for the most 
dominant pairing channel even though the band filling is 
considerably away from half filling and t' is large. 

This result provides an example of singlet pairing ver- 
tex being larger than that of the triplet pairings even for 
an open shell filling, which conversely supports our argu- 
ment in section H*V Al that APf > APd in the quasi-lD 
system may not simply be an artifact of the open shell 
filling, and that /-wave may be competitive against d- 
wave even in the on-site Hubbard model due to effects 
beyond RPA-type theories. 



B. Closed shell filling 

In Fig|51 we present results for a closed shell filling. 
The interaction vertex for spl pairing becomes larger 
than for the open shell filling (FigEJ. The vertices for 
triplet pairings this time take negative values for most 
of the distances. This is in contrast with the result for 
the quasi ID system in Fig[Sl where A/ remains positive 
even for a closed shell filling. This comparison further 
enforces our argument that /-wave pairing is indeed a 
favorable pairing in the quasi- ID system. 



C. Near full filling 

To make sure that we have been on the right track, 
we present here QMC results for a nearly full filled (di- 
lute hole density) case, where the outcome can be safely 
predicted. Namely in this case, spin triplet pairing is ex- 
pected to dominate strongly over singlet since ferromag- 
netic spin fluctuations arise for t' ~ 0.5 and n ~ 2 (or 
equivalently t' ~ —0.5 and n ~ 0) due to the large den- 
sity of states at the Fermi level. This can be considered 
as well established because a number of approaches, in- 
cluding a FLEX studji^ii^, perturbational studies 49 i 50 , a 
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FIG. 7: The QMC result of AP a for a N = 12 x 12 2D system 
with an open shell filling close to n ~ 4/3. (a) focuses only 
on the triplet channels, while (b) shows all the channels con- 
sidered. Parameter values are (N e , n, f') = (190, 1.32, 0.464). 
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FIG. 8: A plot similar to FigEJbut for a N = 14 x 14 system. 
Parameter values are (N e , n, i')=(258, 1.32, 0.406). 



T-matrix approximations^, and a renormalization group 
studjiSS, reach (at least qualitatively) the same conclu- 
sion in this parameter regime. 

Since pairs may be formed at large distances for dilute 
career density, here we calculate the pairing interaction 
vertex also for the singlet sp3 channel shown in the in- 
set of Fig llOf b). For an open shell filling, the interaction 
vertices are in fact positive for all three triplet channels, 
while they are small or negative for the singlet spl and 
sp3 channels as seen in Fig llOf a). By contrast, for a 
closed shell filling, although the interaction vertices for 
tpl and tp3 remains positive, they are much smaller than 
that for the singlet sp3 channel. This result contradicts 
with the established result, again supporting our argu- 
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FIG. 9: A plot similar to FigEJb) but for a closed shell filling. 
Parameter values are (N e , n, t', ty)=(194, 1.35, 0.465, 0.999). 
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merit that the result AP^ > APf for a closed shell filling 
in the quasi ID system (FigJHJ) may not be taken directly. 



FIG. 10: The QMC result of AP a for a 12 x 12 2D system 
near full filling for an open shell (a) and for a closed shell (b). 
Parameter values are (N e , n, t')=(254, 1.76, 0.400). t y = 
0.999 is taken in (b). The sp3 pairing is shown in the right 
inset of (b). 



D. Comparison with FLEX 

To close this section, we present FLEX^Si calculation 
results for the 2D on-site Hubbard model with n = 4/3, 
t' = 0.4 and U = 4, and compare them with the QMC 
results given above. The formulation is as follows. First, 
the renormalized Green's function is obtained within the 
FLEX method, which is kind of a self-consistent RPA. 
Secondly, the pairing interaction mediated by spin fluc- 
tuations is calculated within the RPA form. Finally, the 
FLEX Green's function and the pairing interaction are 
substituted for the linearized Eliashberg equation, which 
is solved using the power method to give the largest eigen- 
value and the corresponding eigenfunction (the gap func- 
tion) . 

The gap functions having the largest eigenvalues are 
shown in Fig llll for (a) the triplet and (b) the singlet 
pairing channels. The singlet gap function has a form 
close to cos(k x ) — cos{k y ), which corresponds to spl pair- 
ing in the QMC study, while the form of the triplet gap 
function is close to sin^a, + k y ), which corresponds to 
tp2. Although the eigenvalues are small (less than 0.3 
for T > 0.005) for both singlet and triplet, the ratio be- 
tween the singlet and triplet eigenvalues turns out to be 
Asingiet / Atripict — 5 at T — 0.005, which is at least quali- 
tatively consistent with the QMC results. 



(a) (b) 




FIG. 11: Contour plots of the gap functions for (a) the triplet 
and (b)the singlet channels obtained by FLEX + Eliashberg 
equation for the 2D Hubbard model with t' = 0.4, n = 4/3, 
(7 = 4, and T(temperature in units of t)=0.005. The thick 
lines represent the nodal lines of the gap function. 



VI. SUMMARY 

To summarize, we have studied the competition among 
various pairing symmetries in the quasi-lD and the 2D 
on-site Hubbard models by calculating the pairing inter- 
action vertices using the ground state QMC technique. 
For the quasi-lD system, /-wave pairing not only domi- 
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nates over p-wave in agreement with the spin-charge fluc- 
tuation theory, but also looks competitive against d-wave 
pairing, /-wave being competitive against d-wave in the 
on-site Hubbard model is rather unexpected from the 
spin-charge fluctuation theory because 2k p charge fluc- 
tuation is not enhanced in the on-site Hubbard model4^ 
If /-wave is competitive against d-wave even in the ab- 
sence of strong charge fluctuations, it is likely that /-wave 
dominates over d-wave in the actual (TMTSF^X, where 
2k F CDW actually coexists with SDW,±2*±2 

For the 2D system, the calculation results (at least for 
the open shell filling) support presence of attractive in- 
teraction in the triplet pairing channel. However, singlet 
d x 2_ y 2 pairing dominates over all the triplet pairings con- 
sidered in the present study at the band filling n ~ 4/3, 



which corresponds to Sr2Ru04, in agreement with the 
FLEX calculation. Although we cannot completely rule 
out the possibility of other (longer ranged) triplet pair- 
ings dominating over singlet, or a possibility of triplet 
dominating over singlet for larger system sizes and/or 
larger values of U, the present results do suggest that it 
is worth considering effects beyond the single band Hub- 
bard model, including the contributions from the a — (3 
bands in order to fully understand the occur- 
rence of triplet superconductivity in Si'2Ru04. 

Numerical calculation has been performed at the fa- 
cilities of the Supercomputer Center, Institute for Solid 
State Physics, University of Tokyo, and at the Computer 
Center, University of Tokyo. 
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